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ABSTRACT 

The nature of dust in distant galaxies is not well understood, and until recently few direct dust measurements 
have been possible. We investigate dust in distant star-forming galaxies using near-infrared grism spectra of the 
3D-HST survey combined with archival multi-wavelength photometry. These data allow us to make a direct 
comparison between dust towards star- forming regions (measured using B aimer decrements) and the integrated 
dust properties (derived by comparing spectral energy distributions [SEDs] with stellar population and dust 
models) for a statistically significant sample of distant galaxies. We select a sample of 163 galaxies between 
1.36 <z< 1.5 with Hg SNR > 5 and measure B aimer decrements from stacked spectra. First, we stack spectra 
in bins of integrated stellar dust attenuation, and find that there is extra dust extinction towards star-forming 
regions (Av,hii is 1.81 times the integrated A v, star ), though slightly lower than found for low-redshift starburst 
galaxies. Next, we stack spectra in bins of specific star formation rate (log sSFR), star formation rate (log SFR), 
and stellar mass (logM*). We find that on average Ay,Hn increases with SFR and mass, but decreases with 
increasing sSFR. The amount of extra extinction also decreases with increasing sSFR and decreasing stellar 
mass. Our results are consistent with the two-phase dust model - in which galaxies contain both a diffuse and 
a stellar birth cloud dust component - as the extra extinction will increase once older stars outside the star- 
forming regions become more dominant. Finally, using our B aimer decrements we derive dust-corrected Hg 
SFRs, and find evidence that SED fitting produces incorrect SFRs if very rapidly declining SFHs are included 
in the explored parameter space. 

Subject headings : dust, extinction — galaxies: evolution — galaxies: high-redshift 


1. INTRODUCTION 

While dust makes up only a very small fraction of the bary- 
onic mass in galaxies (Draine et al. 2007), it can leave a large 
signature on the spectral energy distributions (SEDs). Dust 
extinguishes light in a wavelength-dependent way, which dis- 
torts the intrinsic SED of most galaxies. The dust extinction 
as a function of wavelength must be well understood before 
the intrinsic SED can be recovered from observations. The 
total dust attenuation in a galaxy also depends on the dust ge- 
ometry. Developing a complete model of dust properties and 
distribution in galaxies thus requires observations of both the 
integrated dust properties and the dust towards star-forming 
regions. 
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The properties of dust in low-redshift galaxies have been 
extensively studied. The integrated dust attenuation (A V star , 
which affects the continuum light in a galaxy) has been mea- 
sured with three general methods: (i) the Lir/Luv ratio (also 
known as IRX), which probes dust attenuation using energy 
conservation. This ratio is directly related to the UV contin- 
uum slope P (Meurer et al. 1999), which is used to infer the 
dust content for galaxies for which no IR data is available. 
The other methods are (ii) SED or color analysis using stel- 
lar population models, and (Hi) line of sight analysis using a 
luminous background object (see Conroy 2013 for more dis- 
cussion on this topic). 

Dust extinction towards Hu regions (Av,hii) is most directly 
probed with recombination line ratios, often using the Balmer 
decrement (Hg/H/ 3). The intrinsic line ratio can be calculated 
given reasonable environmental parameters. As dust extinc- 
tion is wavelength dependent, the measured line ratio com- 
pared with the intrinsic ratio combined with an assumed dust 
law yield a measure of the amount of dust attenuation towards 
star-forming regions. This method was used by Brinchmann 
et al. (2004) and Garn & Best (2010) to measure extinction 
towards star- forming regions for SDSS galaxies. 

By comparing the two, Calzetti et al. (2000) find there is 
extra dust extinction towards star-forming regions relative to 
the integrated dust content for local starburst galaxies. Wild 
et al. (2011) expand on this by finding the amount of extra 
extinction increases with the axial ratio of galaxies, implying 
that the dust content of galaxies might have two phases (e.g. 
Calzetti et al. 1994; Chariot & Fall 2000; Granato et al. 2000): 
an optically thick component associated with the star-forming 
regions and a diffuse component distributed throughout the 
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ISM. 

At higher redshifts, the nature of dust is much more poorly 
understood. Most work on dust in high-z galaxies has focused 
on the UV slope (e.g. Reddy et al. 2006; Reddy et al. 2010; 
Wilkins et al. 2011; Bouwens et al. 2012; Finkelstein et al. 
2012; Reddy et al. 2012; Hathi et al. 2013, also see Shapley 
2011 for a comprehensive review), as this is relatively easy to 
observe. However, deviations from the Meurer et al. (1999) 
IR X-/3 relation have been found for galaxy samples with dif- 
ferent properties (e.g. Kong et al. 2004; Johnson et al. 2007; 
Conroy et al. 2010; Gonzalez-Perez et al. 2013), so using the 
UV slope to measure the continuum dust attenuation likely 
requires careful calibration. Indirect methods for measuring 
the attenuation towards Hu regions have also been explored, 
relying on metallicity (e.g. Sobral et al. 2012), IR luminosity 
(e.g. Ibar et al. 2013), and comparison of SFR indicators (e.g. 
Wuyts et al. 2013). 

Direct dust measurements toward Hu regions using B aimer 
decrements are very challenging for z > 0.5 as both Ha 
and H (3 are shifted to the less-accesible near-infrared win- 
dow. Careful survey design and instrument improvements 
have made measurements of the B aimer decrement possible 
for larger samples of intermediate redshift galaxies, e.g. be- 
tween 0.4 < z < 1 (Savaglio et al. 2005), z ~ 0.5 (Ly et al. 
2012), and z ~ 0.8 (Villar et al. 2008; Momcheva et al. 2013). 
However, until recently B aimer decrements have been mea- 
sured for only a small number of more distant galaxies (e.g. 
Teplitz et al. 2000; van Dokkum et al. 2005; Hainline et al. 
2009; Yoshikawa et al. 2010). 

Interestingly, current studies of dust properties in distant 
star-forming galaxies yield contrasting results. There is ev- 
idence for no extra extinction towards star-forming regions 
from Erb et al. (2006b) and Reddy et al. (2010), as extra ex- 
tinction leads to an overestimate of the Ha SFR relative to 
the UV slope SFR. However, other studies find evidence that 
there is extra extinction. Forster Schreiber et al. (2009) com- 
pare measured and predicted L Ha and Ha equivalent widths 
and find the best agreement when Ay 5 hh includes extra extinc- 
tion relative to A v?star . A comparison of overlapping objects 
with Erb et al. (2006b) finds that the previous aperture correc- 
tions might be overestimated, which could mask some extra 
extinction. Yoshikawa et al. (2010) compare Ay 5 star from SED 
fitting and Ay 5 Hii from B aimer decrements for a small sam- 
ple and find the high-z objects are consistent with the local 
universe Calzetti et al. (2000) prescription for extra dust at- 
tenuation. Additionally, Wuyts et al. (2011b; 2013) find the 
best agreement between Ha SFRs and UV+IR SFRs when 
extra extinction is adopted, either the same as the local re- 
lation (Wuyts et al. 2011b) or a slightly lower ratio (Wuyts 
et al. 2013). Mancini et al. (2011) also compare SFRs and 
find the best agreement between Ha SFRs and SED SFRs is 
obtained when the local universe prescription for extra dust 
extinction is adopted. Recently, Kashino et al. (2013) use SFR 
comparisons and the B aimer decrement for a relatively large 
(although targeted) sample and find the amount of extra ex- 
tinction is lower than the local relation. 

These contrasting results may not be surprising, given the 
different and indirect methods and/or small and biased sam- 
ples of most studies. Direct measurements of a statistical sam- 
ple of distant galaxies are required to clarify these dust prop- 
erties. This is now possible, as new IR instruments with mul- 
tiple object spectroscopy capabilities are able to measure the 
B aimer decrement for even larger and more complete sam- 
ples of high-redshift objects. In particular, the Hubble Space 


Telescope's WFC3/G141 grism filter provides slitless spectra, 
allowing for a non-targeted survey of a large number of high- 
redshift galaxies. The HST grisms also avoid atmospheric 
near-IR absorption. 

A number of surveys, including the 3D-HST survey (van 
Dokkum et al. 2011; Brammer et al. 2012), CANDELS 
(Koekemoer et al. 2011), and the WISP survey (Atek et al. 

2010) , have taken advantage of the HST grism capabilities to 
survey high redshift galaxies. Dominguez et al. (2013) were 
the first to use WFC3 grism data to make measurements of 
the Balmer decrement on a large, non-targeted sample. How- 
ever, as their sample was not drawn from regions of the sky 
with existing deep photometric coverage, they were unable to 
examine trends of dust versus integrated galaxy properties. 

We present a statistical study of dust extinction measured 
using the Balmer decrement for a large, non-targeted sam- 
ple of galaxies at z ^ 1.5. Both rest- frame optical spectra and 
deep photometry are available, allowing us to compare extinc- 
tion towards Hu regions with a number of integrated galaxy 
properties. 

Throughout this paper we adopt a ACDM universe with 
D m = 0.3, Da = 0.7, andF/o = 70 kms -1 Mpc -1 . 

2. DATA 

2.1. Observations and catalog 

Our sample is drawn from the 3D-HST survey (Brammer 
et al. 2012), a Hubble Space Telescope (HST) Treasury pro- 
gram adding ACS and WFC3 grism observations to the well- 
covered CANDELS (Koekemoer et al. 2011) fields: AEGIS, 
COSMOS, GOODS-S, and UDS. The 3D-HST data also in- 
clude observations of the GOODS -N field from program GO- 
11600 (PI: B. Weiner). 

The 3D-HST survey provides spatially resolved slitless 
grism spectra with ACS/G800L and WFC3/G141. In this 
work we only use the WFC3/G141 grism spectra, which 
covers 1.1 /im < A < 1.65/im. The raw grism dispersion is 
46.5 pixel -1 , but interlacing during data reduction improves 
the dispersion to ^23 pixel -1 . The G141 grism has a max- 
imum resolution of R ~ 130, corresponding to ~ 110 in the 
middle of our wavelength range. 

The 3D-HST survey makes use of existing deep photomet- 
ric coverage in each of the survey fields, combining the grism 
spectra with the publicly available, multi-wavelength pho- 
tometric catalogs of AEGIS and COSMOS (Whitaker et al. 

2011) , GOODS-N (Kajisawa et al. 2011), GOODS-S (Wuyts 
et al. 2008), and UDS (Williams et al. 2009). 

A modified version of the Eazy code (Brammer et al. 
2008) is used on the combined grism spectra + photometry to 
measure the redshifts, emission line fluxes, and rest-frame U, 
V, J fluxes of the 3D-HST galaxies. Stellar masses, integrated 
dust extinction, SFRs, and sSFRs are determined by fitting 
stellar population synthesis models to the photometric data 
using the FAST code (Kriek et al. 2009). We use a separate 
set of parameters than those used by Brammer et al. (2012), 
for reasons discussed in Section 4.4. We use the Bruzual & 
Chariot (2003) stellar population synthesis models, assuming 
a Chabrier (2003) stellar initial mass function, solar metallic- 
ity, an exponentially declining star formation history with a 
minimum e-folding time of log 10 (r m i n /yr) = 8.5, a minimum 
age of 40 Myr, and an integrated dust extinction Ay between 
0 and 4 assuming the dust attenuation law by Calzetti et al. 
( 2000 ). 

Brammer et al. (2012) provide complete details on the 3D- 


Dust measurements of SF galaxies at z ~ 1.5 in 3D-HST 


3 


o 


DD 

LL 

U3 


cn 

o 


3 


2 


1 


0 


-1 



9.0 9.5 10.0 10.5 11.0 

log 1D (IVL/M 0 ) 


;> 

3 



m 

o 


1.0 - 


0.5 - 


0.0, 


- 0.5 


Fig. 1.- — Sample characteristics relative to all 3D-HST galaxies in the same redshift range. The left panel shows log M* vs. log SFR (from SED fitting) for 
all 3D-HST galaxies at 1.36 < z < 1.5. The black circles indicate the sample selected for direct Balmer decrement measurements, consisting of galaxies with 
a strong H a detection (i.e. S/N > 5). Colors indicate the stellar Ay of the best-fit SED models found using the grism spectra and the photometry, as described 
in Section 2.1. The grey dashed lines show constant values of log sSFR. The right panel presents the rest-frame U-V and V-J colors for all 3D-HST galaxies 
between 1.36 < z < 1.5 (in grey) and our sample (color-coded by stellar Ay bin). We discard any galaxies lying within the quiescent box (using the definition 
from Whitaker et al. 2012a) from our sample, as the emission lines for these galaxies likely originate from AGNs. In both panels, objects within the quiescent 
box of the UVJ diagram are shown as grey squares, while star-forming objects are shown as circles. 


HST survey data reduction and parameter measurement pro- 
cedure. 

2.2. Sample selection 

We select galaxies in the redshift range 1.36 < z < 1.5, for 
which Hct and H/3 are generously covered by the G141 grism. 
In addition we impose a SNR cut of Ha SNR > 5, to mea- 
sure a decent line signal. We make no H/3 SNR cut, to avoid 
biasing our sample against the dustiest galaxies. 

We have a number of additional selection criteria, to ensure 
high quality of the spectra. First, to avoid line misidentifi- 
cation, the photometric and grism redshifts must have good 
agreement, i.e. |z p hot _ %ism| < 0.2. Second, the contamina- 
tion from other sources may not exceed 15% (an issue because 
of the slitless nature of the grism spectra). Third, there must 
be grism coverage of at least 95% (which must include the 
Ha and H (3 lines). Fourth, no more than 50% of the spectrum 
may be flagged as problematic (due to bad pixels or cosmic 
rays) during reduction. 

To study dust attenuation towards star forming regions, we 
do not want AGN to contaminate our emission lines. To reject 
AGN, we exclude any objects that have a detected X-ray lu- 
minosity L x > 10 42 erg s -1 (Mendez et al. 2013; Rosario et al. 
2013; Bauer et al. 2002) by matching against the Chandra 
Deep Field North and South surveys (Alexander et al. 2003; 
Xue et al. 2011) and the XMM-Newton serendipitous survey 
(Watson et al. 2009). Furthermore we reject all objects that 
fall within the Donley et al. (2012) IRAC AGN region, or that 
fall within the quiescent box in the UVJ diagram defined in 
Whitaker et al. (2012a) (as the line emission likely originates 
from an AGN). 

Finally, we visually inspect the grism spectra and photom- 
etry of the preliminary sample to reject problematic objects 
(i.e. wrong line identification). Our final sample includes 163 
galaxies. Figure 1 shows how our sample compares to the full 


galaxy distribution at a similar redshift. The selected galaxies 
all have relatively high SFRs, and lie along the “star-forming 
main sequence” (e.g., Noeske et al. 2007; Daddi et al. 2007; 
Wuyts et al. 201 la; Whitaker et al. 2012b; Nelson et al. 2013). 
The values for our sample fall in the range 0 to 2.3 for Ay star , 
-10.44 to -8.11 for log sSFR, -0.25 to 2.01 fo6r log SFR, and 
8.87 to 11.24 for log 10 M*. 

2.3. Stacking 

The spectra of individual galaxies in this sample are too 
noisy to get a clear measurement of the Balmer decrement, 
the ratio of the flux of Ha to H/3 (F Ha /F H /3). Thus we bin 
galaxies by parameter (SED Ay, sSFR, SFR, stellar mass) and 
stack the spectra before measuring line fluxes. 

We adopt a uniform methodology for stacking spectra 
within a bin. First, we only use the portion of the grism spec- 
tra that falls between 1.13/im and 1.65/im (observed wave- 
length) to avoid noise at the edge of the grism coverage. Then 
the individual spectra are continuum normalized by scaling 
the bi weighted mean value of the flux between 5500 and 
6000A (rest-frame) to unity. The individual spectra are then 
interpolated onto a common rest-frame wavelength grid. 

The normalized spectra are stacked and averaged based on 
the number of spectra contributing at each wavelength. The 
stacked flux errors are simulated using a Monte Carlo error 
estimation, assuming the errors are normally distributed on 
either side of the mean. 

The best-fit stellar population models are sampled over the 
same wavelength regime as the original spectra and used as 
the continua. The continua are stacked in the same way: they 
are continuum normalized and interpolated onto a common 
rest- frame wavelength grid, then averaged together. 

The spectra and photometry have a slightly different slope 
for most galaxies, which results in a tilt between the stacked 
spectra and the stacked continua. To correct for this tilt offset, 
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Fig. 2. — Stacked spectra for bins in stellar Ay (top), log sSFR (second row), log SFR (third row) and log M* (bottom). In each panel, the stacked spectrum 
is plotted in black, with the continuum fit shown in red. The best fit line measurements (from shortest to longest wavelength: H/3, [OIII] AA4959, 5007, blended 
Ha+[NII]AA6548,6584, and [SII] AA6717, 6731) are shown in orange. The continuum fitting is done using the portions of the spectra within the shaded grey 
regions. 


we fit a second-order Legendre polynomial to the continuum 
portions of the stacked spectrum (shown as the grey regions 
in Figure 2). We also fit the best-fit stellar population model, 
while masking the B aimer absorption features. The two poly- 
nomial fits are used to correct the tilt between the stacked 
continuum model and the stacked spectrum. The stacked, tilt- 
corrected continuum is subtracted from the stacked spectrum. 
This removes both the underlying continuum from our spec- 
trum, and corrects for the intrinsic B aimer absorption of stel- 
lar atmospheres. The stacked spectra and tilt-corrected con- 
tinua for the different bins are shown in Figure 2. 

The adopted continuum normalization scheme leads ob- 


jects with higher scaled Ha fluxes (the flux we measure af- 
ter normalizing the spectrum) to have more weight in our line 
stack. To correctly compare the parameters from SED fitting 
(Ay, star? stellar mass, SFR, sSFR) with values calculated from 
the stacked lines, we compute the weighted average of each 
parameter. We use the scaled Ha fluxes as the individual ob- 
jects’ weights. Errors on the average are estimated with boot- 
strapping. 

2.4. Line measurement 

To measure line fluxes, we fit each line in the continuum- 
subtracted stacked spectrum using least- squares minimiza- 
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tion. For our sample redshift range, the grism spec- 
tra have rest-frame coverage of the following spectral 
lines: H/3, [OIII]AA4959,5007, Ha, [NII]AA6548,6584, and 
[SII]AA6717,6731. However, the resolution of the grism data 
is insufficient for us to separate the Ha and [Nil] lines, so we 
measure the blended Ha+[NII]AA6548,6584 line. 

The grism line shapes are not well described by gaussian 
profiles. Because the spatial resolution of the WFC3 detec- 
tor (~0". 13 pixel -1 ) is much greater than the spectral reso- 
lution, the spectral line profiles are dominated by the object 
shapes. Thus, for each object we measure the object shape 
by summing the direct image (a postage stamp extracted from 
the F140W or F160W direct image) over the spatial direction 
(perpendicular to the dispersion direction) and use this as the 
line profile. The composite line profiles are created by flux- 
normalizing the individual profiles, multiplying each profile 
by the object’s scaled Ha flux (described in Section 2.3), and 
finally averaging. This method yields a composite profile with 
the same effective weighting of the objects that results from 
the spectrum stacking method. 

The grism spectra have roughly constant spectral resolu- 
tion, so the composite line profiles described above have the 
correct width for an observed wavelength of 1 .4 /im, or ap- 
proximately Ha in the rest- frame for our redshift range. Thus 
for each line, we scale the line profile width by \n n& / \n a . 

Lines in a spectrum may not have the same profile, possibly 
due to dust or age gradients (e.g. Wuyts et al. 2012). However, 
the line fits obtained while using the same line profile (with 
appropriate width and amplitude scaling) match the data very 
well, suggesting that using a single profile for a stack is a 
reasonable approximation. 

Because of the low spectral resolution of the grism spectra, 
we simultaneously fit the [OIII] doublet and H/3, and simi- 
larly the blended Ha+[NII] line and the [SII] doublet. We fix 
the line ratio between [OIII]A5007 and [OIII]A4959 to 3:1 to 
reduce the number of degrees of freedom in our fit, and we fix 
the redshift of all lines to the value measured for Ha. 

We compute the line fluxes from the best-fit line profile pa- 
rameters. We then perform Monte Carlo simulations to esti- 
mate the flux errors. A simulated spectrum is generated us- 
ing the continuum-subtracted stacked flux and errors. The 
continuum is remeasured and subtracted from the simulated 
spectrum, to include the error from the tilt-corrected contin- 
uum in the fits. The line fluxes are measured on the simulated 
spectrum following the same procedure detailed above. This 
process is repeated 500 times, and the upper and lower 1 a 
errors are calculated from the set of fluxes measured from the 
simulations. 

The best fit line measurements for our stacks are shown in 
Figure 2. 

2.4.1. [Nil] correction 

To measure the B aimer decrement, we need to correct the 
blended Ha+[NII] line for the [Nil] contribution. We use the 
stellar mass versus [Nil] A65 84/Ha relation measured in Erb 
et al. (2006a) for galaxies at z ~ 2, as our sample covers a 
similar range of masses and SFRs, and are close in redshift. 

The stellar masses Erb et al. (2006a) used are based on the 
integrated SFH, not the current stellar mass which we use. For 
the galaxies in our sample, which are all reasonably young, 
the mass from the integrated SFH is about 10% higher than the 
current stellar mass. We use this estimate to scale down the 
masses given by Erb et al. to match our stellar mass definition. 

We interpolate the values Erb et al. report in Table 2 to esti- 


mate the ratio of [Nil] A65 84/Ha with errors for the weighted 
average stellar mass in each bin. We assume an intrinsic line 
ratio of 3:1 between [NII]A6584 and [NII]A6548 to scale the 
ratio to include the second [Nil] line. We use the interpolated 
and scaled ratio to calculate the Ha flux in each stack. 


3. DUST EXTINCTION COMPARED WITH GALAXY PROPERTIES 
3.1. Measuring dust extinction towards star-forming regions 


The Balmer decrement, Ha/H/3, lets us determine the 
amount of dust extinction towards star-forming regions by 
comparing the measured ratio with the expected line ratio 
given the physical conditions of the region. To calculate the 
dust extinction Ay hh from the Balmer decrement, we follow 
the same method as used in Dominguez et al. (2013) and sim- 
ilar to the one used in Momcheva et al. (2013). 

We assume that the Hu region has a temperature T = 10 4 K, 
an electron density of n e = 10 2 cm -3 , and the ions undergo 
case B recombination, so we assume an intrinsic ratio of 
(Ha/H/3)i nt = 2.86 (Osterbrock & Ferland 2006). We assume 
the reddening curve k( A) of Calzetti et al. (2000), so we may 
write 


E(B-V)=1.971og 10 


(Ha/H/3) obs 

2.86 


( 1 ) 


(i.e. see Dominguez et al. (2013)). We combine this with 
the relation Ay = (4.05 db 0.80) x E(B-V ) from Calzetti et al. 
(2000) to calculate the extinction Ay hh from the Ha and H/3 
flux measured for each stack. 


3.2. Integrated stellar Ay 

We first investigate Ay H n in bins of Ay star , to better con- 
strain the currently contested relationship between the two for 
high-redshift galaxies. We choose bins of Ay star to probe the 
full range of integrated stellar dust extinction in our sample, 
from low to medium to high extinction. We define our bins as: 
0 < Ay < 0.6, 0.6 < Ay < 1.2, 1.2 < Ay < 3.0. We stack the 
spectra in these bins and measure Ay hii on the stack using the 
relations given in Section 3.1. The results are shown in Figure 
3 . 

We perform a least-squares ratio fit to the data with a sig- 
nificant (> 2(j ) H/3 detection. The best-fit relation is 

Ay Hu = ^81-043 Ay star , (2) 

indicating that on average Ay hh is 1.81 times higher than 
Ay star in star forming galaxies at z ~ 1.5. This is a slightly 
lower amount of extra extinction than the ratio of 2.27 which 
Calzetti et al. (2000) find for low redshift star-forming galax- 
ies, but they are consistent within the errors. Our results are 
inconsistent with the assumption of no extra dust extinction 
towards star-forming regions (e.g., as in Erb et al. 2006b, for 
agreement between Ha and UV SFRs). 

Our finding is consistent with that of Wuyts et al. (2013), 
who find an average relation between Ay star and A extra = 
(Ay Hu) Hn - Ay star for galaxies at 0.7 < z < L5. Their re- 
lation (shown in Figure 3 by the dotted purple line) is the 
dust extinction required for agreement between Ha SFRs and 
UV+IR SFRs, or if there was no IR detection, SED SFRs. 
Our relation is also consistent with the relation Kashino et al. 
(2013) find by comparing UV and Ha SFR indicators, but 
is not consistent with the relation they find using the Balmer 
decrement. 

We also show the results without correcting the Ha flux for 
[Nil] (open grey circles). These points significantly overes- 
timate the amount of Aynn relative to Ay star , demonstrating 
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Fig. 3. — Ay ;St ar vs. Ay Hu measured from the spectra stacked in bins of 
Ay, star ? shown in the top panel of Fig 2. The black circles show the measured 
Ay, hii using the Balmer decrement (Ha/H/3), while the grey open circles 
show the value of Ay hii if the blended Hct line was not corrected to remove 
the [Nil] flux. The size of the data points corresponds to the fraction of the 
total H a weight in each bin. The Ay, sta r errors shown are the la scatter within 
the bins. The blue dashed line shows ratio of Ay hii to Ay ?star from Calzetti 
et al. (2000). The orange dash-dot-dot line shows the ratio of Ay, hii to Ay, star 
assuming no extra dust extinction towards emission line regions. The black 
line shows the best-fit line to our solid data points, which has a slope of 1.81. 
The fit error is shown with the shaded grey region. This indicates there is 
extra extinction towards emission line regions, but not as much as Calzetti 
et al. (2000) finds necessary for low redshift galaxies. Our best-fit line is 
consistent with the findings of Wuyts et al. (2013), shown with the purple 
dotted line. The green long dash and green dash-dot lines show the relations 
found by Kashino et al. (2013). 

the necessity of correcting for [Nil] when measuring the ex- 
tinction towards star-forming regions using grism data. 

However, it is important to note that there may be sig- 
nificant scatter in the A y, star and Ay ? hh values for individual 
galaxies, so this result only holds on average for a collection 
of galaxies. 

3.3. Specific star formation rate 

In this section we probe the change in dust properties over 
bins of sSFR. We select bins of sSFR to have nearly equal 
numbers of galaxies per bin, or as close as possible given the 
SED sSFR values. The bins we use are defined as: -10.45 < 
log 10 sSFR < -9.31, -9.31 < log 10 sSFR < -8.92, -8.92 < 
log 10 sSFR < -8.10. 

The measured Ay, hii for these stacks are shown in the top 
left panel of Figure 4. The weighted average A V star for the 
sSFR stacks is also shown. The left panels of Figure 4 show 
the difference between A^hh and Ay ?star (AAy, middle) and 
the ratio Ay 5 Hii/Ay ?star (bottom) as a function of sSFR. 

There is a small downward trend in Ay star with increas- 
ing sSFR. There is a much more noticeable decrease in Ay hii 
with increasing sSFR. We also find a decreasing trend in AAy, 
which describes the extra extinction towards star-forming re- 
gions. (We may write A Vj hii = A v?star +A v?e xtra, so AAy = 
^v, extra)- However, it is unlikely to have negative extra dust 
extinction towards star-forming regions. The negative value 


of AAy for the highest sSFR value suggests the Ay star values 
for the objects in that bin may be slightly off. To quantify this 
trend, we perform a least-squares linear fit to Ay extra vs. log 
sSFR, and find the best fit relation is 

Ay, extra = -13.99® - 1.60^;| 2 log^sSFR/yr 1 ). (3) 

The ratio Ay hh /A y star is shown for direct comparison with 
the Calzetti et al. (2000) dust law. On average, the ratio is 
below the value they find, though it is consistent with their 
value for the lowest two sSFR bins. The highest sSFR bin 
is significantly lower than their value. As with AAy, there 
appears to be a decrease of the ratio with increasing sSFR. 
This could be explained by the two-phase dust model, as we 
discuss in Section 4.1. 

3.4. Star formation rate 

Next, we investigate the change in dust properties over bins 
of SFR. We select bins of SFR to have nearly equal num- 
bers of galaxies per bin. The bins we use are defined as: 
-0.26 < log 10 SFR < 0.76, 0.76 < log 10 SFR < 1.07, 1.07 < 
log 10 SFR <2.02. 

The results for stacks in SFR are shown in the center panel 
of Figure 4. Both Ay star and Ay hh are roughly constant for 
the lowest and middle SFR bins, and both increase at the high- 
est SFRs. 

The trend in AAy shows a low amount of extra extinction 
for the lowest and middle SFR bins, but it is consistent with 
zero. The extra extinction for the highest SFR objects is closer 
to ~ 0.7 on average. 

Similar to sSFR, we show the ratio for direct comparison 
with past work. On average, the ratio is below the value found 
in Calzetti et al. (2000). There is no apparent trend in the ratio 
as a function of SFR. 

3.5. Stellar mass 

We also measure dust properties over bins of stellar 
mass. Bins of stellar mass are chosen to have equal num- 
bers of galaxies per bin. Thus we define the bins as: 
8.86 < log 10 M* < 9.82, 9.82 < lo g 10 M* < 10.26, 10.26 < 
log 10 M* < 11.25. 

The results for the stacks in stellar mass are shown in the 
right panel of Figure 4. We see an increase in both Ay star 
and A V hii with increasing stellar mass, with the most massive 
galaxies being the dustiest. 

The plot of AAy shows a roughly constant amount of ex- 
tra extinction for the smallest and middle mass bins, with 
Ay, extra ^0.25. The extra extinction for the highest mass ob- 
jects is much larger on average, with Ay extra ~ 1.5. 

Again, we show the ratio for comparison. On average, the 
ratio is below the value found in Calzetti et al. (2000), al- 
though the smallest and largest mass bin ratios are consistent 
with their value. There is no apparent trend in the ratio as a 
function of stellar mass. 

4. DISCUSSION 
4.1. Physical interpretation 

The observed extra extinction towards emission line re- 
gions, and the decrease in the amount of extra extinction with 
increasing sSFR are consistent with a two-component dust 
model (see Calzetti et al. 1994; Chariot & Fall 2000; Granato 
et al. 2000; Wild et al. 2011). This model assumes there is a 
diffuse (but possibly clumpy) dust component in the ISM that 
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Fig. 4. — Ay measurements for bins in log sSFR, log SFR and log M*. The top row shows Ay star and Ay H ii vs. the binned parameters. The middle row 
shows the difference between emission and stellar Ay vs. the binned parameters. The orange dash-dot-dot line shows the case of no extra extinction towards 
star-forming regions. The bottom row shows the ratio of Ay hh to Av ;Star vs. the binned parameters. The blue dashed and orange dash-dot-dot lines in the third 
row are the ratio between emission Ay and stellar Ay used in Calzetti et al. (2000) and the case of no extra extinction towards star-forming regions, respectively. 
The size of the data points is described in Fig. 3. The errors in log sSFR, log SFR and log M* are the la scatter within the bins. In the center left panel the black 
dotted line shows the best-fit relation of Ay extra = Ay hh - Ay star vs. log sSFR, with the fit errors shown with the shaded grey region. 


affects both the older stellar populations and star-forming re- 
gions, as well as an optically thick dust component associated 
with the short-lived stellar birth clouds that only affects the 
stars within those regions. 

Under the two-phase model, we expect the galaxies with 
the highest sSFRs to have the continuum light dominated by 
young, massive stars. These massive stars would predomi- 
nately still reside in the birth clouds, so both the emission 
lines and the continuum features would be attenuated by both 
the birth cloud and the diffuse dust components, resulting in 
A v ? Hu ~A V star . Galaxies with lower sSFRs would have a 
smaller continuum contribution from massive stars, so more 
of the continuum light would come from stars only attenu- 
ated by the diffuse dust, resulting in A v , H n greater than A v , star . 
These different cases are illustrated in Figure 5. 


Our result of extra extinction towards star-forming regions 
is consistent with studies of local starburst galaxies (e.g. 
Calzetti et al. 2000), with the notable difference that the 
amount of extra extinction they measure is larger than what 
we find. In order to understand this difference, we examine 
the extra dust-corrected equivalent widths of the B aimer lines 
as a proxy for sSFR, and find the Ha equivalent widths for the 
galaxies used by Calzetti et al. (2000) (discussed in Calzetti 
et al. 1994) are much smaller than is measured for our stacks. 
This may imply that the Calzetti et al. sample has lower sS- 
FRs than our sample, so we would expect a larger amount 
of extra extinction for their sample than for ours. However, 
Ha equivalent widths may not be accurate sSFR indicators 
for very dusty systems. In addition, comparisons between the 
two samples are not entirely straight-forward, as the Calzetti 


Price et al. 


High Specific SFR Galaxy Low Specific SFR Galaxy 



Av, stars ~ Av, HI I 


* * * 



Fig. 5. — Illustration of the two-component dust model in galaxies with 
high (left panel) and low (right panel) specific SFRs. The yellow regions 
indicate the diffuse dust component in the ISM. The red regions indicate the 
optically thick dust component associated with the short-lived stellar birth 
clouds. The large blue stars show the young, massive stars which mostly 
are found in the birth clouds. The small red stars show the less massive 
stars (both young and old), which are found both within the birth clouds and 
elsewhere. For galaxies with higher specific SFRs, we expect the continuum 
light to be dominated by the young, massive stars in the birth clouds, so both 
the continuum and emission lines are attenuated by both dust components. 
Galaxies with lower specific SFRs would have a higher contribution to the 
continuum emission from less massive stars, which generally reside outside 
the birth clouds and are only attenuated by the diffuse dust component, while 
the emission lines are still attenuated by both dust components. Thus this 
leads to larger differences between Ay ;Star and Ay hit- 


et al. sample lies above the local star-forming main sequence 
(which may be indicative of mergers) while ours lies along 
the star-forming main sequence. 

A more straight-forward and fair comparison is between our 
findings and those of Wild et al. (201 1), as their sample most 
likely covers the full local star-forming main sequence. At the 
highest sSFRs they find similar amounts of extra extinction 
as Calzetti et al. (2000). They also find increasing amounts 
of extra extinction with decreasing sSFRs, in agreement with 
the trend we observe. Wild et al. find higher amounts of extra 
extinction than we do, but the average sSFRs of their sample 
are lower than for our sample. The two-phase dust model 
naturally explains this difference, based on the dependence 
of extra attenuation on sSFR as discussed above (i.e. lower 
sSFRs lead to larger amounts of extra extinction) for distant 
galaxies and observed in local galaxies by Wild et al. (2011). 

Along the same lines, the two-component dust model 
could explain the discrepancies found between different high- 
redshift studies, if the samples consist of galaxies with differ- 
ent ranges in specific SFR. For example, the Kashino et al. 
(2013) SFR-M* relation implies a higher average sSFR than 
our sample, and they find a lower amount of extra attenua- 
tion. Erb et al. (2006a) find evidence for no extra extinction, 
but the average sSFR is even higher than for the Kashino et al. 
sample. 

Our explanation for the trend between extra Ay and sSFR 
was previously mentioned by Wild et al. (201 1). They suggest 
that the trend may alternatively be caused by a decline in dif- 
fuse dust with decreasing sSFR. However, we observe Ay, star 
increases slightly with decreasing sSFR, so we do not expect 
a decline in diffuse dust with decreasing sSFR. 

In absolute terms, we find that Av,hii increases with mass 
and SFR, and decreases with sSFR. As stellar mass and SFR 
are correlated, the trend of increasing dustiness with SFR and 
mass could share the same cause. As the trend with increasing 
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Fig. 6. — Ay hii vs. stellar mass comparison between this work and past 
studies. The filled data (this work; van Dokkum et al. 2005; Momcheva et al. 
2013; Dominguez et al. 2013; Kashino et al. 2013) indicate direct measure- 
ments of Ay hii using B aimer decrements, while the open data indicate Ay hii 
measured with a different method that is calibrated using B aimer decrements. 
The dashed line gives the median relation derived by Garn & Best (2010) 
using SDSS star-forming galaxies. With the exception of the single object 
from van Dokkum et al. (2005) (combined with information from Kriek et al. 
2006), all other data are the result of stacks (using various stacking schemes, 
both mean and median) or are mean (Momcheva et al. 2013, combined with 
stacking) or median values (Garn & Best 2010; Sobral et al. 2012) of individ- 
ual objects. 


mass is stronger, it is likely that mass is the key property. The 
physical cause may be that more massive objects are able to 
retain a higher fraction of produced metals, or the cause may 
be that lower mass galaxies have larger gas fractions, which 
includes a large amount of less enriched gas. 

4.2. Dust extinction vs. axial ratio 

The two-phase dust model also predicts a dependence of 
dust attenuation properties on the axial ratio. Under this 
model, the amount of dust attenuation from the stellar birth 
clouds is similar in face-on or edge-on systems, while the 
longer path length in edge-on systems results in a larger over- 
all A v, star and a smaller amount of extra attenuation towards 
the star-forming regions. Wild et al. (2011) find evidence for 
the two-phase model from the trends of extinction with axial 
ratio for objects in the local universe. 

The spatially resolved WFC3 images yield excellent axial 
ratio measurements. However, the axial ratio distribution for 
our sample is heavily biased towards face-on systems. The 
more edge-on systems may be dustier, so our selection cri- 
teria likely introduce this bias against edge-on systems. It 
might also be that our sample objects are not disk-like galax- 
ies. There may also be a correlation between the completeness 
fraction and axial ratio. Because of sample incompleteness 
and the small range in axial ratios, we are unable to test the 
two-phase model using the inclinations of the galaxies. 

4.3. Comparison of results for Ay vs. stellar mass 
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Fig. 7. — Comparison between measured Ha SFRs (Kennicutt 1998) and 
the SED SFRs. The two sets of SED SFRs differ only in the choice of the 
minimum star formation e-folding timescale, r: log 10 (r/yr) > 7.0 (open red 
squares) and log 10 (r/yr) >8.5 (closed black circles). (The dashed black line 
shows equal Ha SFRs and SED SFRs.) The excellent agreement between our 
Ha SFRs and the SED SFRs calculated with the higher r m / n provide an inde- 
pendent verification that setting r mzn ~ 300 Myr yields the most reasonable 
SED fits for star-forming galaxies. (See Section 4.4 for a full discussion.) 

A number of past studies have measured B aimer decre- 
ments (and often A v ,hh) versus stellar mass (Kashino et al. 
2013; Dominguez et al. 2013; Momcheva et al. 2013; van 
Dokkum et al. 2005). Other studies have employed differ- 
ent methods of measuring A v ,hh versus stellar mass. These 
include comparing Ljr to Ln a as by Ibar et al. 2013, or cal- 
ibrating [OII]/Ha ratio as an alternate to B aimer decrements 
as by Sobral et al. 2012. We compare our results with their 
findings in Figure 6. 

The collection of results over a range of redshifts makes 
it tempting to speculate that there is no redshift evolution in 
the mass-Ay 5 Hn relation. However, the current data are not 
conclusive. First, there is noticeable scatter in the relation 
as well as incompleteness at the high mass end. Second, the 
measurements of Ay 5 Hii are not necessarily equivalent. 

The samples of Dominguez et al. (2013) and Kashino et al. 
(2013) are similar to ours in redshift and method (stacking) 
for measuring Ay. H ii- Our data show some disagreement at 
high masses (although are consistent within the errors), but 
this may be explained by differences in the galaxy samples or 
analysis methods. Concerning the work by Dominguez et al. 
(2013), we first note that their sample contains fewer objects 
than our sample. Second, we have a broad range of deep pho- 
tometric coverage for our objects, which reduces the errors in 
our masses taken from SED modeling. We also do not com- 
bine measurements between different grisms, to avoid possi- 
ble normalization mismatches affecting the line fluxes. Com- 
pared with Kashino et al. (2013), our sample is non-targeted 
and on average has lower sSFRs. The trend we observe be- 
tween Ay ? hii and sSFR shows that the average Ay 5 hh decreases 
with increasing sSFR, which explains why their values of 
Ay, hii are lower than what we measure for similar masses. 


4.4. Implications for Ha SFR compared with SED SFR 

With our direct dust measurements towards Hii regions, we 
can for the first time determine corrected Ha SFRs using the 
relation by Kennicutt (1998) for a large sample of high red- 
shift galaxies. In this section we use these Ha SFRs (taken 
from stacks in SFR) to test the much more debated SED SFRs. 

The comparison of SFR indicators is shown in Figure 7. We 
compare Ha SFRs measured from bins in SFR and the corre- 
sponding average SED SFRs derived with fitting parameters 
that include an exponentially decaying SFH and a minimum 
e-folding time of log 10 (r m i n /yr) = 8.5 (shown as filled black 
circles). For these SED fitting parameters, we find there is 
excellent agreement between the Ha and SED SFRs. 

This is in agreement with the work of Wuyts et al. (201 lb). 
They compare UV+IR and SED SFRs and find excellent 
agreement between the SFR indicators for log 10 (r m i n /yr) = 

8.5. They also compare the SFR indicators when a shorter 
exponential decay time of log 10 (r m i n /yr) = 7.0 is adopted, but 
in that case they find that the SED SFRs underestimate the 
true SFRs. We repeat our analysis, using the same object 
groupings for stacking and identical SED fitting parameters 
with the substitution of log 10 (r m i n /yr) = 7.0 (shown as open 
red squares). When we allow shorter decay times, we also 
find the SED SFRs are lower, sometimes significantly, than 
the Ha SFRs, in agreement with the findings of Wuyts et al. 
(2011b). 

Reddy et al. (2012) also compare UV+IR SFRs with SED 
SFRs. However, instead of using a longer r mxn for a decreas- 
ing SFH, they find the SED SFRs agree with UV+IR SFRs 
when increasing SFHs are adopted. Though increasing SFHs 
may also give similar results, they are not necessary in order 
to produce SED SFRs that are in agreement with the Ha SFRs 
in this work. 

4.5. AGN contamination 

As discussed in Section 2.2, we use a number of methods 
to reject AGNs from our sample, in order to avoid biasing 
our probe of dust specifically in star-forming regions. How- 
ever, as almost all individual objects do not have sufficient line 
SNRs, and more importantly we do not have separate [Nil] 
and Ha measurements, we cannot use a BPT (Baldwin et al. 
1981) diagram to distinguish whether the line emission origi- 
nates from star formation or AGN. 

After stacking, there is sufficient line SNR to place the 
binned values in the BPT and alternative diagrams, which we 
show in Figure 8. These diagrams allow us to assess if there is 
significant AGN contamination in the stacks. We use the data 
from the stacks in stellar mass in these diagrams. 

In the left panel, we show the traditional BPT diagram: 
[Nil] /Ha vs. [OIII]/H/3. Included are the theoretical limit for 
star-forming (SF) galaxies from Kewley et al. (2001) (dashed 
line) and the more conservative empirical division between 
SF galaxies and AGN by Kauffmann et al. (2003) (dashed- 
dot line), both of which are derived for galaxies at z ~ 0. The 
SDSS DR7 galaxies (Abazajian et al. 2009) are divided into 
three categories based on these dividing lines: SF galaxies 
(blue contours), composite (green contours), and AGN (red 
contours). These same group definitions are used in the alter- 
nate BPT diagrams (middle and right of Figure 8). Because 
we do not measure [Nil] directly (instead inferring a value 
as described in Section 2.4.1), we cannot use the traditional 
BPT diagram as a post-analysis check. However, we still can 
see if our inferred [NII]/Ha ratio is sane given the observed 
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Fig. 8. — Left: Traditional BPT (Baldwin et al. 1981) diagram used to identify the ionizing mechanism for optical emission lines. The dashed line is the 
theoretical limit between star-forming (SF) galaxies and AGN, as derived in Kewley et al. (2001), and the dashed-dot line is the empirical division derived in 
Kauffmann et al. (2003) for galaxies in the SDSS. The color contours represent SDSS galaxies, divided into three regions according to these dividing lines. The 
darker, thicker lines enclose 68% of the population, while the ligher, thinner lines enclose 95%. The grey open points represent our stacked spectra (in bins 
of Af*), with the [NII]/Ha measurement inferred from the average stellar masses (i.e. not directly measured from the data.) Middle: alternative BPT diagram, 
using the blended lines observed with the grism. Color contours correspond to the SDSS galaxies in the left panel. The black circles show the measurements 
from our stacked spectra. This panel illustrates that the combination of blended emission lines does not enable us to discriminate between SF galaxies and AGN. 
Right: the Mass-Excitation (MEx) diagnostic from Juneau et al. (201 1), where the solid lines are the empirical divisions between SF galaxies and AGN, valid for 
z < 1 galaxies. Our stacked spectra (black circles) lie a bit above the division between SF galaxies and AGN, but recent studies have shown that high redshift SF 
galaxies tend to be offset up and to the right from lower redshift SF galaxies in this diagram (see Section 4.5). 


[OIII]/H/3 ratio, which we plot as open grey circles to empha- 
size that these are not real measurements. We note that the 
data do lie within the conservative star forming galaxy region, 
so our inferred [NII]/Ha ratios are consistent with values for 
SF galaxies. However, in agreement with other observational 
(e.g. Shapley et al. 2005; Liu et al. 2008) and theoretical stud- 
ies (Kewley et al. 2013) of distant galaxies, the data points are 
moved slightly up and right in the BPT diagram compared to 
the local star- forming sequence. 

An alternate BPT diagram, used in Dominguez et al. 
(2013), is [SII]/Ha vs. [OIII]/H/3. However, Ha is not a 
directly measured quantity, so using the deblended Ha val- 
ues makes it impossible to detangle the assumptions of the 
[Nil] correction with the possible presence of AGN. Instead, 
we make an alternate diagram with only directly measured 
quantities: [SII]/(Ha+[NII]) vs. [OIII]/H /?, which we show 
in the middle panel of Figure 8. There is significant overlap 
between the different categories of SDSS galaxies, so we are 
unable to discriminate between SF galaxies and AGN using 
combinations of blended emission lines. 

A final alternate to the BPT diagram, the Mass-Excitation 
(MEx) diagram (Juneau et al. 2011), is useful as we directly 
measure the necessary lines ([OIII] and H/3) and measure stel- 
lar mass with SED fitting, which we show in the right panel 
of Figure 8. The solid black lines are their empirical divi- 
sions between SF, composite, and AGN (valid up to z ~ 1). 
We include the SDSS categories (using stellar masses from 
the MPA/JHU value-added catalogs; Kauffmann et al. 2003) 
to demonstrate that while the separation is not as clean as the 
traditional BPT diagram, it does separate SF galaxies from 
AGN in the local universe. Our data lie above and to the right 
of the SF region as defined by Juneau et al. (201 1). However, 
recent work at high-z (Newman et al. 2013) has shown that 
the low redshift empirical divisions in the MEx diagram in- 
correctly classify high-z SF galaxies as AGN. They use the 
shift in the mass-metallicity relation between z ~ 0 and 2 to 
derive a shift left (to lower mass) to bring the high-z values 
into agreement with the z ^ 0 MEx diagram (see Fig. 15 of 
Newman et al. 2013). With such a shift, our data shifts se- 


curely into the SF region of the diagram. This demonstrates 
that there is likely not much AGN contamination in our data 
sample. 

4.6. [Nil] contamination 

One of the largest sources of uncertainty in this analysis 
is our correction for [Nil] in our measurement of Ha from 
the blended Ha+[NII] line. Figure 3 demonstrates how much 
our results change when we do not correct for [Nil]. We use 
the relation between [NII]/Ha and stellar mass by Erb et al. 
(2006a) as a way of detangling the blended lines. We in- 
fer the [Nil] contribution from the weighted average mass of 
each stack. Because the two samples have similar masses and 
SFRs, and are close in redshift, the average sample proper- 
ties should be fairly similar. Still, the use of this relation may 
introduce bias into our results. 

If the gas phase metallicities of the Erb et al. (2006a) sam- 
ple are on average lower than for our sample, for example 
due to the fact that their sample is at slightly higher redshift 
(z ~ 2), we may be underestimating the [Nil] contamination. 
Also, it is possible there is some AGN contamination at the 
highest masses (Kriek et al. 2007), which would also result 
in an underestimate of the [Nil] fraction. Finally, the differ- 
ence in stellar mass definition (Erb et al. use the integrated 
SFH, while we use the current stellar mass), even after our 
attempt to correct for the difference in definition, may result 
in discrepant masses and thus biases in the [Nil] correction. 

In most cases, our measured Ha flux would be higher than 
the true value, leading to an overestimate of A v ,hh- Also, 
the integrated stellar Ay and mass are correlated. Thus if the 
[NII]/Ha versus mass relation is incorrect, there may be a sys- 
tematic offset of A y Hu at high mass or A V star . 

4.7. Incompleteness and other systematic uncertainties 

One of the strengths of our analysis is we draw a sample 
from a non-targeted grism survey, with sample cuts designed 
to avoid bias as much as possible. However, bias and incom- 
pleteness most likely still affect our sample. 
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The dustiest star- forming galaxies have very attenuated Ha 
fluxes. Our Ha SNR selection cut, designed to avoid adding 
noise to our analysis, introduces bias against galaxies with 
large Av,hii- Given the galaxy distribution in Figure 1, the in- 
completeness of the most dusty galaxies will result in missing 
the most massive star-forming galaxies, and possibly some of 
the highest SFR galaxies. 

Our continuum normalization scheme also introduces bias 
into our analysis. We adopt a normalization scheme to im- 
prove the signal of our stack, but the cost is that some objects 
have much higher scaled Ha fluxes, and thus they dominate 
our stacks. This biases our results towards those objects with 
higher Ha equivalent widths in a given bin. Because the line 
measurements are biased based on the Ha flux, we take the 
weighted average of the continuum values within a bin to en- 
sure a fair comparison. However, if we instead used a non- 
weighted average, the A v , star vs. A v ,hh relationship does not 
change much. 

Finally, the Ay star derived from SED fitting may suffer from 
severe systematic uncertainties. For example, assuming a 
fixed attenuation law, implemented as a simple screen leads to 
biases in the integrated dust attenuation measurements (e.g., 
Kriek & Conroy 2013). Degeneracies with other modeling pa- 
rameters may also lead to biases. Including rest- frame mid-IR 
data in future work will ensure more accurate values of Ay star . 


5. SUMMARY 

In this paper, we investigate dust attenuation in distant star- 
forming galaxies using data from the 3D-HST survey. We 
measure both the dust towards Hu regions, using B aimer 
decrements, and the integrated dust properties, using SED 
modeling. We find that there is extra extinction towards star- 
forming regions. On average the total attenuation of these 
regions, Aynn, is 1.81 times the integrated dust attenuation, 
Ay star- 

However, the amount of extra attenuation is not the same for 
all galaxies. We find that the amount of extra attenuation de- 
creases with increasing sSFR, in agreement with the results by 
Wild et al. (2011) for low-redshift galaxies. The extra extinc- 
tion also increases with increasing SFR and stellar mass. Our 
findings are consistent with the two-phase dust model, which 
assumes there is a diffuse dust component in the ISM and 
an optically thick dust component associated with the short- 
lived stellar birth clouds. For galaxies with high sSFR, the 
stellar light will be dominated by continuum emission from 
the younger stellar population in the birth clouds, resulting in 
similar attenuation toward the line and continuum emission. 
For more evolved galaxies, much of the stellar light will only 
be attenuated by the diffuse ISM, leading to larger discrepan- 
cies between the two dust measures. 
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Similar to previous studies (Forster Schreiber et al. 2009; 
Yoshikawa et al. 2010; Wuyts et al. 2011b, 2013; Kashino 
et al. 2013), we find less extra extinction in distant galaxies 
than is found in the local universe (Calzetti et al. 2000). This 
effect can also be explained by the two-phase model, as lower 
redshift objects tend to have lower sSFRs than higher redshift 
galaxies (e.g. Noeske et al. 2007; Whitaker et al. 2012b; Fu- 
magalli et al. 2012). 

We also find that both Ay H n and Ay star increase with in- 
creasing SFR and stellar mass, and decreasing specific SFRs. 
This illustrates that more massive galaxies with higher SFRs 
are also more dusty in general. This is generally consistent 
with previous studies of Ay hh and stellar mass over a wide 
range of redshifts. We also observe there to be little redshift 
evolution in the Ay H n-A/* relation, although errors and in- 
completeness makes it impossible to make a definite claim. 

We use our corrected Ha SFRs to test the accuracy of SFRs 
derived from SED fitting. We find excellent agreement be- 
tween the SFR indicators if short SFH decay times are not 
allowed and the constraint log 10 (r/yr) > 8.5 is used. This 
independently confirms the results of past studies comparing 
UV+IR and SED SFRs (Wuyts et al. 2011b) or Ha and SED 
SFRs (Wuyts et al. 2013). 

We note that our results are slightly impacted by incom- 
pleteness and systematic uncertainties. First, we employ SNR 
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could impact our measurements. Finally, grism spectra have 
low spectral resolution, so we deblend the Ha and [Nil] emis- 
sion lines using a metallicity-derived [Nil] correction scheme, 
which could introduce systematic bias to our results as well. 
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tions by a number of new multi-object near infrared spectro- 
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allowing for investigation of dust to higher Ay limits. 
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